Dissipation in graphene and nanotube resonators. 
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Different damping mechanisms in graphene nanoresonators are studied: charges in the substrate, 
ohmic losses in the substrate and the graphene sheet, breaking and healing of surface bonds (Velcro 
effect), two level systems, attachment losses, and thermoelastic losses. We find that, for realistic 
structures and contrary to semiconductor resonators, dissipation is dominated by ohmic losses in the 
graphene layer and metallic gate. An extension of this study to carbon nanotube-based resonators 
is presented. 

PACS numbers: 73.23.-b, 03.65.Yz, 62.25.+g, 85.85. +j 
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I. INTRODUCTION. 

Nano-electro-mechanical devicesii^i^ (NEMS) have at- 
tracted a great deal of attention, as they are a problem of 
fundamental interest, and also because of their potential 
applications. 

Recently, NEMs made from graphene sheets have been 
studied^, following work on NEMs based on carbon 
nanotubes^i^. These devices show unique characteristics, 
as graphene sheet stacks have a high elastic modulus and 
very small total mass thanks to the low number of atomic 
planes (sometimes just even one) composing the bridge 
or cantilever, allowing for higher resonating frequencies 
than other materials of similar dimensions and increased 
potential sensitivity. 

Graphene itself has attracted a great deal of attention^, 
because of its unique features. Graphene samples can be 
made one carbon layer thick, and doped by an external 
electric field. The lattice dynamics of these thin samples 
have not been studied in detail yet. Two dimensional sys- 
tems have, in addition to acoustic modes, transverse flex- 
ural modes^, which are the ones explored in the experi- 
ments considered in this paper. The quadratic dispersion 
of these modes lead to a constant density of states at low 
energy. Experimental observations show that free stand- 
ing graphene is not flat^ i^°'^^ , but corrugated. These rip- 
ples imply the existence of flexural deformations, and can 
lead to charge inhomogeneities in single layer graphene at 
low dopingsi^ii^. Most graphene samples stand on Si02 
substratesiii, and the interaction between the graphene 
layers and the substrate is not well known. The experi- 
ments discussed in this paper can provide information on 
this issue. 

The potential sensitivity of a resonator-based detector 
may be in practice strongly reduced by dissipative pro- 
cesses affecting the vibrational mode used for detection, 
due to the associated widening of the resonance, which 
masks the frequency shifts used to determine the pres- 
ence of external species adsorbed or close to the detector. 
Hence, it is of fundamental importance to gain knowledge 
about those damping mechanisms, to establish their rela- 



tive importance and dependence on resonator parameters 
(dimensions, elastic constants, temperature, etc), which 
may help to optimize performances and determine where 
should efforts be put, not only to use them as detec- 
tors, but also as tools for the study of more fundamental 
questions like the quantum to classical crossover with in- 
creasing system sizes^""'^^. 

In common resonating structures made of semiconduc- 
tors the prevailing dissipative mechanism with decreasing 
size and temperature are surface-related losses: the pres- 
ence of the imperfect surface, with its roughness, struc- 
tural defects, impurities and dangling bonds, can be mod- 
eled by a distribution of effective two-level systems which 
couple to the vibrational eigenmodes of the devic o^^i^"^'^^ . 
Many other processes contribute to a lower extent to the 
damping of vibrations in these devices. Some of them 
are common to all experimental setups, like attachment 
losses 'iSi^S or thermoelastic damping2ii22iiS^. Others de- 
pend on the actuation scheme: For example, in the mag- 
netomotive actuation scheme^!, a layer of metal is de- 
posited on top of the vibrating semiconducting structure 
to control its motion with the Lorentz force actuating 
on the electron current that passes through the top layer 
in presence of an applied magnetic field. This metallic 
layer increases dissipation in two different ways: i) In- 
creasing the local temperature due to electron-phonon 
interactions, thus "feeding" other mechanisms whose ef- 
fect tends to grow with T, ii) Absorbing energy through 
the excitation of electron-hole pairs in the metallic layer 
due to the presence of fixed charges in the substrate sup- 
porting the oscillating structure, which create a potential 
on the electrons moving within the mobile structure that 
is time-dependent, as perceived by the latter. 

This last mechanism has been ignored until now in 
the literature, perhaps due to the small amount of fixed 
charges in typical substrates like single-crystal Si or 
GaAs. But in the case of graphene or carbon-nanotube 
based resonators it must be considered, as it plays a 
much more significant role. This is due to the fact that 
graphene is conducting, and in some actuation setups, 
the control over the graphene layer's motion is through 
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FIG. 1: (Color online) Sketch of the system considered in 
the text. V represents the Coulomb interaction between the 
charged graphene layer and the Si gate. 



the establishment of a capacitive couphng between two 
charged layers, namely the oscillating graphene and a 
doped Si backgate. The number of carriers in both can be 
controlled by an external gatei^. The coupling between 
the charges of both layers, apart from enabling the con- 
trol of the resonator's motion, causes energy losses which 
will dominate at high temperatures. There will be also 
fixed charges in the supporting structure, mainly in the 
Si02 layer located between the graphene and the doped 
Si backgate, absorbing energy too from the resonators 
motion. In this paper we give will give full account of 
these processes. 

In the following, we start analyzing different dissipa- 
tive processes which may be present in devices based 
on graphene. Our calculations should give reasonable 
order of magnitude estimates of the strengths of these 
mechanisms. In Sections II and III we model the ab- 
sorption of mechanical energy due to the charge present 
in the oscillating graphene sheet, the Si02 substrate and 
the Si backgate. In section IV we discuss the role as 
attenuation source played in these peculiar devices by 
the breaking and healing of bonds gluing the graphene 
sheet to the Si02 substrate, a possibility also missing in 
current semiconducting resonators. Whereas for the lat- 
ter prevail surface-related losses, we show in Section V 
that for graphene resonators this friction mechanism is 
highly suppressed thanks to their high degree of crys- 
tallinity. For completeness, we apply previous results of 
other works to estimate the effect of two more friction 
sources present in all setups, namely attachment losses 
and thermoelastic losses, in Section VI. Once these mech- 
anisms have been studied, an extension of the results to 
carbon nanotube resonators is presented in Section VII. 

To make numerical estimates we focus on the devices 
studied ini. The parameters which characterize the av- 
erage oscillator studied there are given in Table [I]. A 
sketch of the system is shown in Fig. [T] . 



II. STATIC CHARGES AT THE AMORPHOUS 
SIO2 SUBSTRATE. 

The graphene sheet can couple electrostatically to 
static charges. These charges give rise to a time depen- 



System properties 


Dimensions 




Thickness t 


10 • 10~® m 


Width w 


10"^ m 


Length L 


lO"'* m 


Height above substrate d 


300 ■ 10"^ m 


Frequency fo 


100 MHz 


Amplitude A 


0.5 nm 


Carrier density pc 


lO'^'cm"^ 


Properties of graphite 


Mass density p^vf 


2200 kg/m^ 


Elastic constants 




E 


10^2 Pa 


V 


0.16 


Debye temperature Qd 


~ 570 A" 


Specific heat Cp 


700 J / Kg. K 


Thermal conductivity k 


390 W / m . K 



TABLE I: Parameters used in the calculations presented in 
the main text, adapted to the systems studied in^^. Bulk data 
taken fronj2^. 



dent potential acting on the electrons of the vibrating 
graphene. The energy is dissipated by creating electron- 
hole excitations in the graphene layer. Static charges 
have been proposed as a source of scattering by the car- 
riers in the graphene^^i^. 

The time-dependent component of the unscreened po- 
tential induced by a charge separated by a distance d in 
the vertical direction from the graphene layer, acting on 
an electron at position r in the graphene layer is given, 
approximately, by: 



y(r,t) 



(|f|2+d2)3/2 



(1) 



where A is the amplitude of the flexural mode, and wq 
its frequency. The dissipation depends on the screening 
by the graphene layer—. A calculation of the damping is 
given in Appendix [X] We find that, for a single graphene 
layer, a single charge gives a contribution to the inverse 
quality factor of: 



1 2h 

Mvt 



kpd > 1 
kpd < 1 



(2) 



where M is the mass of the oscillating sheet, fcp = TTy^pci 
and PC is the density of carriers in the graphene sheet. 
Typical values of this quantity are in the range pc ^ 
10i2cni-2, so that kpd ~ 10^ - 10^ > 1. Eq.® can be 
generalized to a graphene sheet with N layers: 

1 2?i 



(3) 



The suppression with the number of layers is due to the 
increased screening in this system. 

The total contribution to the inverse quality factor is 
obtained by multiplying eqs.® or ([3]) by the total num- 
ber of charges iVch. An upper bound to the density of 
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local charges, deduced from some models for the electric 
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conductivity of graphene- ■ ' is Pch 
ing the parameters in Table [l], we find A'ch 
~ 10"^"'^ at low temperatures. 
This mechanism leads to ohmic dissipation, as the en- 
ergy is dissipated into electron-hole pairs in the metallic 
graphene layer. Hence, the temperature dependence of 
this mechanism is given hy Q~^{T) ^ Q^^{0)x{kT/huJo) 
, and Q-^ - 10"^ at 300 K. 



X ReV^",, {r',u;^)x Imx" [? - f , c^q] } (6) 

The static screening properties, lim^_^o R'ex"(q, 0), of 
the graphene layer and the gate are determined by their 
electronic compressibilities, v'-' and i/*^ respectively. We 
will assume that the distance between the graphene and 
the gate is much larger than the electronic elastic mean 
free path in either material, so that their polarizability 
is well approximated by: 



III. OHMIC LOSSES AT THE GRAPHENE 
SHEET AND THE METALLIC GATE. 

The electrons in the vibrating graphene layer induce a 
time dependent potential on the metallic gate which is 
sometimes part of the experimental setup. The energy 
is transferred to electron-hole pairs created at the gate 
or at the graphene layer. These processes contribute to 
the energy loss and decoherence of electrons in metallic 
conductors near gates22i^. 

The coupling between charge fluctuations in the two 
metallic systems is due to long range electrostatic inter- 
actions. The corresponding hamiltonian is 



H = 1 
2 



(z,r,t)p^(z,r, t) + 



1 



1 Et^w r/920N 2 

2 12 Ka?/ 



.(0,r',i)pG(f)} 



ic 2pMtw 

where the indices G and C stand for the gate and 
graphene layer, respectively. pM is the mass density of 
the graphene sheet, and t,w,E its thickness, width and 
Young modulus, whereas (/)(r, t) represents the vibrating 
amplitude field of bending modes and 11 = dL/d4) is its 
conjugate momentum (L is the Lagrangian). The self- 
consistent screened potentials Vscr{z, r, i), Wscr(0, r, t) are 
calculated as a function of the bare potentials vq{z, r, i), 
vo{Q,Y,t) in Appendix IB] 

As in the case of eq.([T|), the time-dependent part of the 
bare potentials couples the electronic degrees of freedom 
and the mechanical ones through the charge p(r) and 
amplitude of the vibrational mode, A^, and would give 
rise to a term in the quantized hamiltonian of the form 



Hint oc p{v)A^ cx [bl + b^) 



k,k' 



h.c] (5) 



where and p(r) have been expressed in terms of cre- 
ation and annihilation operators of phonons (q, lu^) and 
electrons of a 2D Fermi gas, respectively. 

But a realistic model requires taking into account the 
screening of the potential associated to these charge fluc- 
tuations. In terms of the screened potentials, the induced 
broadening of the mode (q, w^) of the graphene layer can 
be written, using Fermi's Golden Rule, as^S 



a=G,C 



x"(q,^) 



(7) 



where = UpP is the diffusion constant, and is the 
elastic mean free path. The two dimensional conductivity 
is = k^l"'. 

We assume the gate to be quasi two dimensional. This 
approximation is justified when the distance between the 
gate and the graphene layer is much larger than the width 
of the gate. In this situation, the broadening of the mode, 
eq.®, can be expressed as 



T{uo^)^ J d^k\vscr{d,Ku^)\^lmx'^ + Wscr{0,KtOii)\^lmx^ 

(8) 

The screened potentials for a graphene layer oscillat- 
ing in an eigenmode (q, of amplitude A^, have in 
a first approximation only one momentum component. 



are (see Appendix [B|) 



Vscr{d, q,LU^) 



,.(q, cj^)(5(k — q), and these components 



PC^q 



|q| 



Wscr(0,q,t^q) = 



(9) 

where q — |q| and pc is the charge density in the 
graphene layer. The results for T{uj^) and Q~^(wq) can 
be formulated in terms of the total charge in the graphene 
layer, Qc = J d^vpc ^ Lxw x pc- 

In the limit of short separation between the layers, d <C 
L, which is the situation present in current experimental 
setups, one has 



4^2 



I 



G.2 1 



(10) 



The limit I?|q|2 ^ cj for the imaginary part of the suscep- 
tibility of a dirty metal, Imx(q, tf) « wi^/Djqp, has been 
used. The first term in the summation describes losses 
at the graphene sheet, and the second at the gate. The 
associated inverse quality factor, according to eq. (|A6|) . is 
given by 
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To make numerical estimates, we use the parameters in 
Table[I], with i^'^{E) = E/2TTh^vl, vp « 10^ m/s for 
a single layer of graphene, and v'~^ [E) = {N'-f)/vp for a 
stack of N layers^. Carriers in graphene stacks have 
large mobilities'^, and we take D'-'v'-^ « 10'^. Typi- 
cal charge densities for the graphene layer are pc ~ 
lO^^cm"^, leading to a total charge Qc ^ 10'*. For 
these parameters, the contribution of the graphene sheet 
is ~ 10^^. The relative contribution from the gate 
depends on the distance to the graphene sheet. For a 
Si layer with D'^u'^ k, 10^ and at short distances, the 
contribution to the damping from the gate is of the same 
order as that of the graphene sheet. 

Damping is associated to the creation of e-h pairs in 
a metal, which implies that this mechanism is ohmic. 
The inverse quality factor should increase linearly with 
temperature, leading to Q^^ ^ 10^^ at 300 K. 



IV. BREAKING AND HEALING OF SURFACE 
BONDS: VELCRO EFFECT. 



In the fabrication process of the device, the graphene 
flake is deposited on the Si02 substrate, and becomes 
linked to it through hydrogen bonds created by the 
silanol groups (SiOH) present at the substrate's surface. 
When the flake is set into motion, some of this bonds may 
repeatedly break and heal (the Velcro effcct'^^), causing 
dissipation of the energy stored in the vibration. Nu- 
merical estimates are difhcult to make, but nevertheless 
two qualitative arguments showing that its role in the 
damping is probably negligible can be presented: 

i) This mechanism is expected to be temperature inde- 
pendent, in contrast with the strong decrease of friction 
observed as temperature is lowered^. 

ii) The elastic energy stored in a typical graphene oscil- 
lator of lateral dimensions w ~ l/im is about lOeV, when 
the amplitude is ~ Inm. This means about ~ 10~^eV per 
nm^. On the other hand, the energy per hydrogen bond is 
about 10~*eV, and typical radical densities at Si02 sur- 
faces are^ ~ lnm~^. Hence the elastic energy available 
on average for each hydrogen bond is much less than the 
energy stored in the bond. Only rare fluctuations, where 
a significant amount of energy is concentrated in a small 
area will be able to break bonds, and to induce energy 
dissipation. Note, however, that this argument ceases to 
be valid for very large amplitudes ^ 30nm. For higher 
amplitudes, this mechanism can induce significant losses. 



V. DISSIPATION DUE TO TWO-LEVEL 
SYSTEMS. 

This is the typical mechanism for the damping of 
sound waves in insulating amorphous materials^ii^i^. 
An atom or a few atoms can have two nearly degenerate 
configurations. A vibration modifies the energy differ- 
ence between these situations. This mechanism leads to 



the damping of acoustic phonons in amorphous Si02^. 
It is also expected to dominate friction in many NEMs^l. 
We expect the graphene sheet to show a high degree of 
crystallinity, and we will only consider two-level systems 
(TLSs) in the rest of the structure. 

The TLSs can only dissipate energy if they are coupled 
to the vibrating graphene sheet. A possible mechanism 
is the existence of charge impurities associated to these 
defects (fluctuating charges), which are electrostatically 
coupled to the conducting electrons in the graphene. 

We expect this mechanism to be less effective in the 
device considered here than in NEMs made of semicon- 
ducting materials, as now the TLSs reside in the Si02 
substrate, not in the vibrating structure. The coupling, 
arising from long range forces, will be in comparison ac- 
cordingly suppressed, by a factor of order (a/d)", where 
a is a length comparable to the interatomic separation, 
and n describes the decay of the coupling [n — \ for the 
Coulomb potential between charged systems). 

The temperature dependence of the contribution of 
TLSs to is determined by the density of states of 
the modes coupled to the TLSs and the distribution of 
TLSs in terms of their paramenters (tunneling amplitude 
Ag and bias Ag)J^.The hamiltonian describing the cou- 
pling of the effective TLS's and the oscillating graphene 
sheet is given by^^ 

A^ 

k k 

(12) 

where e = ^ (Ag)^ + (Ag)^, is the coupling con- 
stant, which will be strongly suppressed in these de- 
vices as compared to attenuation of acoustic waves in 
amorphous materials, 7 ~ leV x (a/c?)", foj^ represent 
the phonon creation operators associated to the differ- 
ent vibrational modes of a sheet, and X^k'^k^^k ~^ ^^k-^ 
represents the coupling to the strain tensor Uik- There 
are two types, compression modes (longitudinal waves) 
and bending modes. The damping is due to the ini- 
tial transfer of energy from the vibrational mode studied 
by the experimentalists to the TLSs, which in a second 
step transfer this energy to the rest of the modes. The 
properties of the spin-boson model, eq. , are fully de- 
termined by the power-law s of the spectral function^, 
= J2k h^k^o/e\^ S{uj - ujk) ^ aw^~^u;^ where ujk 
is the frequency of mode k, a is an adimensional con- 
stant and LOco is the upper cutoff of the phonon bath. 
For this system, compression modes gives rise to a super- 
ohmic, s — 2, bath, while the bending modes constitute 
an ohmic bath, s = 1, and thus will prevail as a source of 
dissipation at low temperatures^. We will therefore re- 
strict our analysis to the dissipation caused by the ohmic 
component of the vibrational spectrum. 
Applying the method ini^ to the 2D bending modes of 
the graphene sheet, one arrives at J{i-u) — aco, with 



a : 



AS\2pl//(l + ^)3/2(l-^^)i/2 



(13) 
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Here v is the Poisson ratio of graphene. Choosing fairly 
symmetrical TLSs, Ag/e ^ 1, for the parameters in table 
m a ~ 10~^ X [a/dY"^, very small. In ref.'^^ an expres- 
sion is given for the inverse quality factor of a vibration 
damped by TLSs in amorphous insulators, 



EkT 



de du ^ C{uj,T) 

(14) 

where u = A^/e, emax ~ 5 K, and {u\/l — u^) ^ comes 
from the probability density of TLS's in an amorphous 
solid, like SiOj. Q-^iuj,T) is a function of C(cj,T), the 
Fourier transform of the correlation function C{t,T) — 
{<7z{t)az{Q))T- For biased TLSs and a ^ 1 an extensive 
analysis of C{lo,T) is performed in^^, where several ex- 
pressions are provided in different limits. Using them, 
the estimate for (5~^(a;,T) follows: 



1 Q-\lu,T) « ^fkT, kT < e,na. 

(15) 

In the range of temperatures of current experiments 
(5K<T<300K), the dependence of dissipation with T is 
weak, and Q^^ - 10"^ x (a/d)'^ - 10^^^. The main 
uncertainty of the calculation has been the use of the 
TLSs' distribution assumed for amorphous solids^, but 
due to the small value of a a weak dissipation is expected 
also with a modified distribution. Thus the conclusion is 
that the relative importance of TLSs damping is much 
smaller for graphene than for other NEMs devices^i^i^. 



B. Thermoelastic effects. 

When the phonon mean free path of the acoustic 
phonons is shorter than the wavelength of the mode un- 
der study, the acoustic phonons can be considered a dissi- 
pative environment coupled to the mode by anharmonic 
terms in the ionic potentia l^^i^^'^'^ . These anharmonic 
effects are described by the expansion coefficient, a, and 
the thermal conductivity, k. We follow the analysis in^. 
For a rectangular beam vibrating at a frequency lu the 
inverse quality factor is 



Qz\T) = 



Ea^T LOTz 



(17) 



where E is the Young Modulus, Cp is the specific heat 
at a constant pressure, and tz is the thermal relaxation 
time associated with the mode, which in the case of a 
flexural vibration is given by tz = t^Cp/{'K^K). This 
estimate assumes that the graphene sheet is weakly de- 
formed, and that the typical relaxation time is associated 
to the diffusion of phonons over distances comparable to 
the thickness of the sheet. 

Although better approximations are available in the 
literature^^, eq. ([TT]) is enough for an estimate of the order 
of magnitude of Q'^ . Using the parameters from table[I], 
for t = 10 nm and / ~ 100 MHz, we find that utz < 1, 
and 



Q'z\T = 300K) 



Ea^TujV 



5 • 10 



-7 



(18) 



VI. OTHER FRICTION MECHANISMS 
A. Attachment losses. 

The energy is transferred from the resonator mode to 
acoustic modes at the contacts and beyon d^^i^° . 

The main expressions needed are given in2£. When 
t, and d is much smaller than the wavelength of the 
radiated elastic waves in the Si02 substrate, the contri- 
bution to the inverse quality factor is given by 




pg,ii;g(i~(^o)2) 



(16) 



where the superscript O applies to the silicon oxide, and 
stands for Poisson's ratio. The range of values of the 
quality factor varies from Q^^ « 5 • 10~^ for a graphene 
monolayer, to « 5 • 10^^ for a stack with 30 lay- 
ers and t — lOnm. These quantities probably overesti- 
mate the attachment losses, as they do not include the 
impedance at the Si02-graphene interface. 

This damping process due to energy irradiated away 
from the resonator should not depend on temperature. 



VII. EXTENSION TO NANOTUBE 
OSCILLATORS. 

The analysis presented here can be extended, in a 
straightforward way, to systems where the oscillating 
part is a nanotube. 

We expect in these devices a larger impedance between 
the modes of the nanotube and those of the substrate, so 
that attachment losses will be suppressed with respect to 
the estimate presented here for graphene. 

The damping mechanisms which require long range 
forces between the moving charges in the nanotube and 
degrees of freedom of the substrate (fluctuating and static 
charges) will not be significantly changed. A nanotube 
of length L at distance d from the substrate will interact 
with a substrate area of order {L + d) x d. A similar 
estimate for a graphene sheet of length L and width w 
gives an area ^ (L + d) x (w + d). As L ~ w d ~ 1/im, 
the two areas are comparable. 

On the other hand, ohmic losses induced in the nan- 
otube will be reduced with respect to the two dimensional 
graphene sheet, as the number of carriers is lower in the 
nanotube. 

Finally, we expect a longer phonon mean free path in 
the nanotube, which implies that thermoelastic effects 
will be reduced . 
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Q""(r = 300A') 


Temperature 
dependence 


Charges in the Si02 


10" ' - 10"" 


r 


Charges in graphene sheet 


10"^ 


T 


aiiLi iiieLdinc ^dLc 






Velcro effect 


Absent 




Two-level systems 




A + BT 


Attachment losses 


IQ-'* - 10"^ 




Thermoelastic losses 


10-^ 


T 



TABLE II: Contribution of the mechanisms considered in the 
main text to the inverse quality factor Q~^{T) of the systems 
studied 



VIII. CONCLUSIONS. 

We have considered six possible dissipation mecha- 
nisms which may lead to damping in a graphene meso- 
scopic oscillator. The main results are summarized in Ta- 
ble [11]. We expect that the calculations give the correct 
order of magnitude and dependence on external param- 
eters. 

We find that at high temperatures the leading damping 
mechanism is the ohmic losses in the metallic gate and 
the graphene sheet. This effect depends quadratically 
with the total charge at the graphene sheet, which can 
be controlled by the gate voltage. 

At low temperatures attachment losses limit the qual- 
ity of the vibration. If the resonator is strongly driven, 
a new damping mechanism may come into play, the 
Velcro effect, which may limit substantially the qual- 
ity factor as compared with the slightly driven case. 
The high crystallinity of the resonators eliminates the 
main source of dissipation in semiconducting resonators, 
namely surface-related effective TLSs coupled to the local 
strain field. 

These conclusions apply with only slight modifications 
to carbon nanotube-based resonators. 
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where xo is the susceptibility of the graphene layer. At 
low energy and momenta its value tends to the compress- 
ibility of the electrons in the layer: 



¥^ Re[xo(|q|,t^)] 

|q|— >0,Ci;— >0 




(A3) 



where N is the number of layers and 7 is the interlayer 
hopping element. For a stack with N layers, we have 
used the model with one interlayer hopping elemen"b^, 
which gives rise to 2N low energy bands, most of which 
show a quadratic dispersion. 

Using Fermi's golden rule, we finally find for width of 
the graphene mode (t;(q, aj)=u(— q, w)): 



Fph- / d^q|t;(q)pImxo(q,Wo) 



(A4) 



where: 



Imxo(q, t^o) 



^ll¥ 



N = 1 



v§.\<i\^/P ^ 



(A5) 



where, for 7^ 1, p is the total carrier density. This last 
expressions are valid for lengths bigger than the mean 
free path, I ^ Imfp- 

The energy absorbed per cycle of oscillation and 
unit volume will be AE = (2Tr/u!o)fiiUQTpY^/twL ~ 
2iThTph/twL, and the inverse quality factor Q~^{uJo) will 
correspond to 



Qph('^o) 



1 AE 



ph 



1 



2hr 



ph 



twL 



(A6) 



where E^ is the elastic energy stored in the vibration, 
M is the total mass of the resonator, and A the am- 
plitude of vibration. Substituting eas. (jAllA2IA3IA5l) in 
ea. (jA4p . and inserting (jA4p in eq. (jA6p . one arrives at 
cqs.([2]) and ^ for the dissipation due to a single charge 
in the substrate. The analysis presented here does not 
consider additional screening due to the presence of a 
metallic gate. In that case, one needs to add to the po- 
tential from a static charge, eq.(IT]) in the main text, a 
contribution from the image charge induced by the gate. 
This effect will reduce the coupling between the graphene 
layer and charges in the vicinity of the gate. 



APPENDIX A: COUPLING TO FIXED 
CHARGES IN THE Si02 SUBSTRATE 

The Fourier transform of the potential in eq.([T])is: 

T/(q, u) = 2Tre^ Ae-'^'^S{uj - ujq) (Al) 
This potential is screened by the polarizability of the 
graphene layer~, so that has to be replaced by : 
^2 1^1 



l + e2/|q|Re[xo(|q|,^)] 



Re[xo(|q|,c^)] 

(A2) 



APPENDIX B: SCREENING OF THE 
POTENTIALS AT THE GRAPHENE SHEET 
AND SI GATE 



The equations for the selfconsistent potentials 
c,(2:,f — f,w) as a function of the bare potentials 
f — f ,Lu) are given by 
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Vscrid,r -r',uj) = Vq {d,r - r' ,lu) + Vq {d,r - r' ,lu) + / dvi / dr2Vcouiid,r - ri,uj)x'^' (vi - r2,uj)vscr{d,r2 - r' ,uj) + 

Jc Jc 

+ drs driVcouiid,r-r3,Lo)x^{r3-r4,uj)vscr{0,ri-r',Lo) 
Jg Jg 

Vscr{0,r - r' ,uj) = VQ{0,r - r' ,oj) + Vq {0,r - r' ,u)) + / dri / dr2i;coui(0, r - Fi, w)x'^(fi - f2, w)wscr(0, f2 - r', w) + 

Jg Jg 

+ drs dr4t;coui(0,r - r3,a;)x*^(r3 - r4,a;)wscr((^,r4 - fjw) , (Bl) 
Jc Jc 



r 



where for example in the first equation Vq {d, r ~ r' , u) graphene layer due to the presence of charges in the Si 

represents the bare potential experienced by a point plane, i^coui is the two-dimensional bare Coulomb po- 

charge e in the graphene layer due to the presence of tential. These equations simplify considerably in the q 

charges in that same layer, while r — r',u}) is the space: 

bare potential experienced by a point charge e in the 

I 



Vscr {d, q, w) = {d, q, a;)e«'* + {d, q, lo) + VgX (q, i^)vscr {d, q, w) + I'^x (q, '^)Vscr (0, q, w) 
Vscr{0, q, w) = v§{0, q, w) + v§{d, q, a;)e«'' + Vqe'i'^x''^ {q, ui)vscr{d, q, w) + w,x*^(q, w)t;scr(0, q, w) 

I 



(B2) 



where Vg = 27re^/|q| is the Fourier transform of of (d, q, w) and Vq {0,q^,Lo). Now we can calculate 
the Coulomb potential in two dimensions, and where Vscr{d,q,,i^) and Wscr(0, q, w) in terms of the rest of the 
Vq{0, q, oj) and VQ{d, q, w) have been expressed in terms variables, 



Vscr{d) \ _ 
Vscr{0) 



1 - VqX^ -VqC-l'^x'^ 

-Vge-f^X^ 1 - VgX^ 



1 e 



v§{d) 



(B3) 



r 



The dependence on q and u) has been omitted for the long wavelenght limit ygx"" ,Vqx'^ ^ 1, the last equation 
sake of clarity. Now, if we are interested only in the simplifies to 



^scr 
Vscr 



id) 
(0) 



v„\x e 



C „-qd 



^G^qd 



C 



Vq [-X^ + x"") Vq (x^e-?'^ - x^e-i") 



r 



«o^(0) 
v^{d) 



(B4) 



a. Values of Vq {0,q,Lo) and Vq {d,^,Lj) 

Now we will calculate the parts of these terms which 
will give rise to a coupling to the vibration. When the 



graphene layer is set into motion with a bending mode of 
wavevector q and amplitude A^, the potential of a point 
charge e in the Si plane due to the charge in the graphene 
layer, v^{0,r,t), is 
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27re^po 



(f-r')2 + (d + A^|e*(^P'-"9<))2 



c 



rf2 



c 



((?-?') 



d2 



3/2 



where in the second hne an expansion for small A^^ has Similarly, the potential of a point charge in the oscil- 

been performed. The Fourier transform for 7^ is lating graphene sheet due to the charge in the Si plane 



(B6) 



is 



v'iid, r,t)^U df , « /(f) + ^e^poA^e^W--.^*) (B7) 



2 



G 

r 



leading to the same expression as eq. (|B6|) but without Substituting (jB6IB8[) in ea. (|B4p . one obtains eq.©. 
the factor e~'^'^ 

i;o^(d,k,c.')=«o''(0,k,^')e'' (B8) 
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